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.  Abstract  (Continued) 

parabolic  models  of  the  ionosphere,  and  the  corresponding  real  height  values 
were  then  obtained  by  the  present  method.  The  results  obtained  showed  that 
the  calculated  values  of  real  height  were  in  very  good  agreement  with  the 
exact  values  of  real  height. 


To  indicate  the  improvement  obtainable  with  the  introduction  of  additional 
parameters,  a  two-parameter  function,  which  gives  a  somewhat  more  accu¬ 
rate  representation  of  the  dispersion  curve  than  does  the  one-parameter 
function,  was  obtained  for  one  case.  However,  the  amount  of  computational 
labor  required  is  considerably  greater  than  that  required  for  the  present 
method. 


In  view  of  its  simplicity  and  accuracy,  it  is  felt  that  the  one-parameter 
method  is  a  useful  analytic  tool  for  determining  real  height  from  ionosonde 
data. 
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A  Simplified  Procedure  for  Direct  Determination 
of  Real  Height  From  Virtual  Height  Data 


1.  INTRODUCTION 

in  a  previous  report1  a  method  was  presented  for  directly  calculating  real 
height  from  virtual  height  data  obtained  in  the  presence  of  a  magnetic  field.  The 
method  was  based  upon  the  representation  of  the  non-linear  portion  of  the  disper¬ 
sion  relation  between  index  of  refraction  and  electron  density  by  a  simple  power 
law,  resulting  in  an  Abel  integral  equation  for  the  relation  between  real  height  and 
virtual  height.  The  equation  may  then  be  solved  explicitly  for  the  real  height  as  a 
function  of  frequency. 

In  making  numerical  applications,  it  was  found  that  the  use  of  different  analyti¬ 
cal  forms  for  the  two  regions  of  the  dispersion  curve  resulted  in  fairly  tedious 
calculations.  In  addition,  taking  account  of  the  variation  of  the  several  parameters 
with  frequency  (in  particular  the  position  separating  the  linear  and  non-linear  por¬ 
tions  of  the  curve!  considerably  enlarged  the  extent  of  the  calculations.  In  seeking 
an  alternative  method,  it  was  found  that  the  entire  dispersion  curve  could  be  rep¬ 
resented  with  fairly  good  accuracy  by  a  single  power  law,  which  greatly  simplified 
and  reduced  the  required  numerical  work. 


(Deceived  for  publication  5  March  1981) 

1.  Klein,  M.M.  (1979)  A  Method  for  Direction  Determination  of  Real  Height  From 


Two  parabolic  models  were  chosen  for  the  investigation:  In  Model  1  the  elec¬ 
tron  density  increases  to  arbitrarily  large  values;  in  Model  2  the  electron  density 
reaches  a  maximum  value.  Model  1  was  chosen  because  of  the  ease  of  making  the 
required  calculations;  Model  2  was  utilized  because  it  has  been  applied  frequently 
in  studies  of  the  ionosphere. 


2.  DETERMINATION  OK  KKAL  HEIGHT 


The  dispersion  curve  is  written  in  the  form 


p“  (1  -  X)1' 


(1) 


where  p  ts  the  index  of  refraction  and  X  is  the  ratio  f^/f~,  where  f  is  the  plasma 
frequency  and  f  is  probing  frequency.  The  exponent  a  is  a  parameter  to  be  deter¬ 
mined  by  fitting  the  curve  given  by  Eq.  (1)  to  the  exact  dispersion  curve.  The 
virtual  height  h'(f),  which  is  a  function  of  input  frequency  f,  is  related  to  the  real 
height  h  by  the  integral 


h'(f)  = 


I 

/ 


g'  dh 


(2) 


where  h  is  the  height  of  reflection,  that  is,  when  g  -  0,  and  p'  is  the  group  index 
efraction  and  given  by 


nr  f) 


Ci) 


We  may  now  write  Eq.  (2)  in  the  form 


d 

ar 


pf  dh  =  h'(f) 


(4) 


in  which,  since  p  vanishes  at  the  upper  limit  h  we  have  taken  the  derivative  out¬ 
side  the  integral  sign. 


Integration  of  Eq.  (4)  yields 


7 


or 


f 

o  v 

where  Fv  or  f  is  the  maximum  value  of  F  or  f  corresponding  to  the  reflection 
height. 

To  determine  the  value  of  a,  the  approximate  dispersion  curve  was  compared 
with  the  exact  curve  for  several  values  of  a.  In  general,  the  approximate  curve 
was  above  the  exact  curve  for  small  values  of  X  and  below  at  the  larger  values. 

A  value  of  a  was  chosen  that  was  estimated  to  give  the  least  overall  error.  A  plot 
of  a  as  a  function  of  frequency  is  given  in  Figure  1  for  two  values  of  0,  where  0  is 
the  angle  between  the  ray  and  the  magnetic  field.  A  comparison  of  the  exact  and 
approximate  dispersion  curves  is  given  in  Figures  2  and  3  for  the  same  two  angles 
The  exact  curve  was  computed  for  several  values  of  the  gyromagnetic  frequency 
Y  =  fj_j/f  at  fjj  =  1.  5  MHz.  It  is  noted  that  the  curves  are  generally  in  good  agree¬ 
ment  at  an  angle  of  30°,  while  at  10°  the  agreement  is  not  as  good,  particularly  at 


Figure  1.  Variation  of  Parameter  a  With  Frequency 


ti 


t 


Figure  3.  Comparison  of  Exact  and  Approximate  Dispersion  Curves  for 
Several  Frequencies;  9  =  30°,  fj  -1.5  MHz 


the  lower  frequencies.  However,  since  the  error  changes  sign  as  X  goes  to  larger 
values,  it  is  anticipated  that  the  calculated  values  of  real  height  for  10°  should  not 
differ  much  from  those  for  30°. 
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The  cltH'iron  density  N  is  ft'vcn 
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The  virtual  height  is,  from  Eq.  (4),  given  by 


( 13) 


h'(f) 


(if 


(f  I) 


where 


l  -  J  (1  -  X)‘>  ;2  dh  . 

O 

Using  Eq.  (12),  we  can  write  I  in  the  form 

if 


(14) 


i  -  2  4  f  u  -  *>tt/2  dx 


(U>) 


where  b  =  f/fp, 
tron  density  Np. 


wMl  ,  ,h,  plasm.  frequency  corresponding  to  (he 

Integration  of  Eq.  (O)  yields 


1 

2  +  a/2 


(16) 


The  virtual  height  h'(f)  is  then,  from  Eq.  (13) 


*»'<n  "  ^ 


(17) 
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A  hu  ll  Yields 


h 1  ( f ) 


r  L?l_  llii 
1  2  df 

c 


(1 


2  -  „ 

- r, - 

„  2)”(2 


9 


2)- 


(18) 


An  examination  of  the  second  term  shows  that,  because  of  the  small  value  of 
(see  Figure  1),  it  is  typically  only  a  few  percent  of  the  first  term.  We  shall  there¬ 
fore,  for  simplicity,  neglect  this  term. 

I  sing  Kq.  (7),  the  quantity  ^  ^  -■  is  given  by 


cl  1 1  '(f) 
elf 


<  If) 


whirl)  yields 


dll’(f)  .  Ji+u 

df  ^  ' 


(20) 


A 


10 _ 1 _ 1 

lt.i'2  2  t-  o' 2  A  2 
f  c 
P 


(21) 


in  which,  in  view  of  the  small  value  of  o  -  1,  we  have*  ignored  variation  in  f1'"1  and 
have  again  neglected  terms  in  . 

The  real  height,  Kq.  (10),  now  takes  the  form 


h 


sin  tt,-, .  2 
~ 72 


rH" 


<fj  -  rf"2 


df 


(22) 


which  can  be  expressed  in  terms  of  a  new  variable,  u 


f/f  , 
v 


a  s 


h 


sin  tfi>/2 
*r,/2 


3+o 


(1 


du 


(23) 


The  integral  in  Kq.  (23)  can  be  expressed  in  terms  of  the  Gamma  function  and 
Kq.  (23)  written  as 


h 


sin  tt <y / 2 

*o/2 


r(l  -  t»/2)r(2  +  e/2) 


(24) 


12 


4.  ANALYSIS  FOR  MODEL  2 


Here  the  electron  density  is  specified  by 


1  - 


?  = 


(25) 


(26) 


where  hm  is  the  height  corresponding  to  the  maximum  electron  density  and  a  is  a 
constant  determined  by  the  value  of  N  at  some  particular  height.  The  virtual  height 
is  then  given  by 


h'(f)  =  4r  <W) 


(27) 


J  = 


q  -  x)^2 

(1  -  b2  X)l/2 


dX  . 


(28) 


The  integral  J  cannot  be  evaluated  exactly  for  arbitrary  b.  For  values  of  b  less 
than  one  a  power  series  solution  may  be  obtained,  while  for  b  close  to  one  a 

combined  numerical  and  analytical  technique  may  be  utilized. 

2  -1/2 

For  b  less  than  one,  the  quantity  (1  -  b  X)  '  can  be  expanded  into  a  power 
series  and  the  integral  J  in  Eq,  (28)  takes  the  form 


J 


a  b 
2 


T  .  1 .2  T  3  .  • 4  5  .  6  3  5,8  . 

J°  ?  Jl+8b  J2  lTb  J3+mb  J4 


(29) 


where 


J"'/ 


(1  -X) 


a/2 


Xn  dX 


n  =  0,  1,  2, 


The  integral  <Jn  can  be  evaluated  by  successive  integration  by  parts,  or  by  express¬ 
ing  in  terms  of  the  Beta  function 
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R(m,  n) 


J  x"'-1 


(1  -X)  dX 


r(m)r(n) 
r  (m  +  n) 


The  integral  J  can  then  bo  written  as 
”  n 


Jn  "  B(n  +  1  *  J  +  X) 

r(n  +  i)r(|  +  lj 
r(n  4-  2  +  a/2) 


which  can  be  reduced  to 


n  1  +  a/2  12  +  a/2  3  +  a/ 2  n  +  1  +  a/2 


1+0/2 


n  1,  2,  3, 


The  virtual  height  h 1  ( f)  is  then  given  by 


a  (oh2T  ,5,4  21  ,6.  45  ,8^  385,  10. 

h’(n  =  -  l3b  Jo  2  b  J  1  T  b  J2  T6  b  J3  128  b  J4  + 


Values  of  the  virtual  height  can  be  obtained  with  good  accuracy  and  with  only  a  Few 
terms  for  b  S  0.  8. 

The  real  height  is,  using  Eqs.  (4)  and  (10),  now  given  by 


sin  tta/2 


rv  f°_1 

f  J  ,.2  :f2.a/2  h'<f>  df 

O  '  V 


which  may  be  expressed  in  the  form 


h  I  3JoKo+!JlKl+fJ2K2  +  llJ3K3  +  li  J4K4  + 
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To  check  the  accuracy  of  the  g  term,  the  value  J  given  by  Eq.  (41)  was  com¬ 
pared  with  that  given  by  Eq.  (29)  for  b  in  the  range  0.  7  to  0.  85;  agreement  to  with¬ 
in  a  few  percent  was  obtained.  The  quantity  is,  using  Eq.  (19),  now  given 

by 


dHJSfl  ,  sjb  ^ 


(43) 


and  the  real  height  h  is  now  calculated  from 


o 


-1  b  in  (\±*j  -j -JfJ? 


(44) 


which  can  be  written  in  the  form 


“  ■  ¥  >■, 


In 


1  +  b  u 
v 


du 


1  -  bv  u/  (1  -  u2)o/2 


(45) 


where 

»  =  f/fv 


The  integral  in  Eq.  (40)  can  be  evaluated  by  expansion  of  the  logarithmic  term 
as  a  power  series  in  b^u.  However,  this  series  expansion  offers  no  advantage  over 
the  series  expansion  in  powers  of  b,  Eq.  (29),  a  large  number  of  terms  being  re¬ 
quired  for  bv  close  to  unity  in  either  case.  Since  the  integral  can  be  evaluated 
exactly  in  the  neighborhood  of  the  singularity  (bv  =  1),  we  have  utilized  numerical 
integration  over  about  95  percent  of  the  range,  and  evaluated  the  remaining  5  per¬ 
cent  analytically.  We  write  the  integral  L  in  Eq.  (45)  as 


(46) 


(47) 
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'•9 


/  u‘"'  (^v) 


du 


(1  -  u2)°  2 


(48) 


and  Uj,  iho  point  separating  I.j  and  l.^,  >-s  about  0.95.  The  integral  1.  was  eval¬ 
uated  numerically  by  Simpson's  rule.  To  evaluate  J.„,  we  change  the  variable  to 
e  1  -  u,  €j  1  -  Uj,  and  obtain 


—™2  (Mj  -  -M2  -  M3) 


(49) 


where 


l-o/2 

fi 

M,  -  In  (1  +  b  )  ~r - s 

l  v  t  -  r_ 


(50) 


bv  *1'*  ~ 
A|  _ • _ .  _ i - 

*U2  1  +  b  2  -  a  '2 


(51) 


In  (1  -  b  +  b  e) 

iM3  /  - rr£ — —  de 


/ 


(52) 


For  b  1,  3,1,.  becomes 

v  ’  3 


M, 


f  In  e 

J  7 72 


de 


1-0/2 

fjL _ 

1  -  o/2  l1"  C1  1  -  a/2 


In  e 


1 


(53) 


For  b  *  1,  1U  cannot  be  expressed  in  terms  of  elementary  functions,  but  can 
v  '*  2 

be  expressed  in  terms  of  the  digamma  function.  Integration  by  parts  yields  for 


M. 


3 


..  ln(I  '  bv+  Vl>  l-o/2  b  M 

1  -  a/2  fl  '  1  -  o,/2  M4 


1 


<53  d 


2.  Abramowitz,  A.,  and  Stegun,  I.,  Editors  (1964)  Handbook  of  Mathematical 
Functions,  National  Bureau  of  Standards,  AMS  55,  Washington,  D.  C,, 
Chap,  fi. 
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( 


2 -u  2  l  .  _  ,  ... 

ei _  r  t1'1-'-  fj _  r  - 

-  b  /  b  c ,  I  1  -  b  /  /  b  ,  e . 

°  1  +r^ir  1  + 1  - '  r^ir 


2-i.  2  1 


l-i. '2 


M 


4  1 


1  +  1 

-  t 

u .  ± 

l  1 

-K) 

dt  (54) 


where  t  e/e^.  Wo  may  now  approximate  by 


2-o  12 


M 


4  ‘  1  -  b 


1 

/ 


tl-«/2 

1  +  t 


dt  +  1 


Vi  ^ 

a. 


t> 


2  d' 


(55) 


where  the  integrals  in  Eq.  (55)  are  given  by 


■  &)] 

■b 

k  >  0 


(56) 


i 

/ 


tk_l  .  1  k  -  1 

5  dt  -  7T - ;; — 

(t  +  1)  2  2 


m  -  -  ® 


(57) 


k  >  0 


and  i/<(z)  is  the  digamma  function 

ii(z)  =  ^  In  r(z)  .  (58) 

5.  A  METHOD  FOR  IMPROVING  ACCURACY  OF  APPROXIMATE 
DISPERSION  CURVE 

As  indicated  in  Section  2,  the  accuracy  of  the  calculated  values  of  real  height 
depends  to  a  large  extent  upon  how  well  the  dispersion  curve  can  be  represented  by 
a  single  parameter  function  such  as  the  form  given  by  Eq.  (1).  The  results  ob¬ 
tained,  as  shown  by  Figure  2,  indicated  good  accuracy  at  30°,  the  accuracy  being 
only  fair  in  the  neighborhood  of  10°.  An  increase  in  accuracy  can  be  achieved  by 
introducing  additional  parameters,  the  significant  practical  restriction  being  the 


3.  Grobner,  W.,  and  Hofreiter,  N.  (1965)  Integraltafel.  Zweiter  Teil.  Bestimmte 
Integale,  Springer-Verlag,  Wien,  New  York,  p,  176. 
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amount  of  computational  labor  introduced.  Because  of  the  non-linearity  of  the 
method,  additional  terms  cannot  be  introduced  in  a  straightforward  manner,  l  or 
example,  in  the  two-parameter  function 


■in 


(1  -  X) 


"1 


+  (1  -  X) 


rt2  2 


(59) 


it  is  tedious  to  obtain  t>j  and  o9  but,  more  significantly,  a  non-integral  power  of 
frequency  is  introduced  into  the  electron  density  term  in  Eq.  (4)  that  render  s  pro¬ 
hibitively  difficult  the  inversion  of  the  integral  equation  for'  real  height.  However, 
a  form  that  does  lead  to  a  practical  solution  is  the  two-parameter  function 


p  =  (1  -  X)a/2  (1  -  aX) 


(90) 


where  a  and  a  are  parameters  to  be  determined.  Integration  of  Eq.  (4)  now  takes 
the  form 


(F-  5)a/2  u  d?  -  a-^F  -  ?)?/2  =  H'(F) 


(61) 


which  can  be  written  as 


(1  -  a)  J  (F  -?)o/2udf  +  -|  J 


F  /it 

(F-'-^  u  d?  =  H'(F) 


(62) 


The  Laplace  transform  of  Eq.  (62)  is 


ds  --  H(s) 


(63) 


w 


hich  can  be 

dv(s) 

ds 


written  as 


v(s)  _ 
s 


1  1  H(s) 

1  -  a  1  +  aj 2  s 


(64) 


where 


v(s)  = 


ds 


(65) 
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The  solution  of  Eq.  ((14)  is 


v(s) 


s_1'b  H(s)  ds 


(Ml) 


where 


b  -  j— - —  (1  +  a  1 2) 

1  -  a 

1  1 

c  "  1  -  a  T(1  +  aJ2)  • 
Differentiation  of  Eq.  (G6)  yields 


dv(s)  _ 
ds 


£.  Il(s)  +  cbsb_1  J 


s'3'1  H(s)  ds 


((17) 


which,  in  terms  of  u(s),  takes  the  form 


u(s)  =  — — 7 w  H(s)  -  cb 

_  1  -of  £ 


<y  OC 

A-*”2 1 


H(s) 

1+b 


ds 


((18) 


The  inverse  transform  of  Eq.  ( fl 8)  is  then 
F 


dh 


-V  2 

c  d~  f  H’(F)  el>  iL 

d?  r(l  -  n/2)  df2  J  (F  .  r)ol2  '  ni  .  b  .  2)  df 

O  V 


2 1 


G(F) 


(F  -  F) 


b-f-o  /  2 


dF 


(H9) 


where 


G(F) 


l- 

_1 _  J,  f 

r  ( 1  -  b)  F  J 


(F  -  t)  H’(t)  dt 


(70) 
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Integration  of  Eq.  ((>9)  and  carrying  through  the  remaining  differentiation  yields 


W2)  / 


1  dH'(F) 


J  (F  -  F) 
o  v 


aj2  dF 


T(1  -  b  -  <r/2) 


(Fv  -  F) 


b+o/2  dF 


.  dG(F)  .  , 

where  — jp-1  is  given  by 


dG(F) 

dF  '  *  T(1 


^7  / 


(F  -  t)b  H'(t)  dt  + 


r 

_ b_ _  f 

r(i  -  b)  J 


(F  -  t)*5'1  H'(t)  dt 


An  approximate  solution,  somewhat  simpler  in  form  than  Eq.  (71),  can  be  ob¬ 
tained  if  a  is  small  compared  to  unity.  Subsequent  numerical  calculations  show  that 
this  is  indeed  the  case.  Neglecting  the  term  in  a  in  Eq.  (63)  and  solving  for  u(s)  in 
terms  of  H(s)  yields 


-  s1+"/2  1  - 

u(s)  =  (1  _  -)  r(l  +  a/2)  H(s)  ’ 


which  enables  us  to  now  write  Eq.  (63)  in  the  form 


77(a)  -  _ 1 _ I  s2  a(l  +  a/2)  s2 

(  (1  -  a)  r(l  +  o/2)  g l-o’/ 2  H(S)  1  -  a  gl-c 


?  00  _ 

r  H(s> 

o/2  J  ~~s 


Taking  the  inverse  transform  of  Eq.  (74)  yields 


dh  sin  Ttorl 2 


d?  (1  -  a)  na/2 


(F  -  fW2 

V 


4  <-  11(f>-;2  dF-ajft! 4  [ 

d?2  J  (F  -  F)a/2  (1  ‘  a)  dF2  J 

o  v 
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where 


To  show  the  improvement  in  the  approximation  to  the  dispersion  curve,  the 
case  0  -  10°,  f  =  5  was  studied  numerically.  The  parameters  a  and  a  were  deter¬ 
mined  by  fitting  the  approximate  curve,  Eq.  (60),  to  the  exact  dispersion  curve  at 
the  points  and  X 2,  setting 

=  (1  -  -  aXj) 

U2  =  (1  -  X2)q/2(1  -  aX2) 


where 

X1 

-  0.  3 

,  /it  =  0.  872 

X2 

0.7 

,  m  2  ■  0.  87  1 

Elimination  of  a  from  these  equations  yields 


In  n  In  ( 1  -  Xj) 

In  q,  in  (1  -  \~) 


(7't) 


)  O 


where 


A  solution  to  Kq.  (70)  was  obtained  graphically,  yielding  a  0.21;  the  value  of 
a  was  then  determined  as  0.40.  A  check  of  the  approximate  dispersion  curve 
against  the  exact  curve.  Figure  3,  now  shows  excellent  agreement  up  to  x  0.  3; 
for  X  >  0.  8  the  improvement  is  only  fair.  It  is  likely  that  a  better  approximation 
for  X  can  be  obtained  by  choosing  a  larger  value  for  X.,,  say  X7  -  0.  8  or  0.  U. 


6.  KKSl  I.TS  \NI>  DISCI  SSION 

The  calculated  values  for  virtual  height  as  a  function  of  frequenc  y  are  presented 
in  Figure  4,  and  those  for  the  real  height  are  shown  in  Figure  5.  The  virtual  and 
real  height  have  been  normalized  to  the  maximum  real  height.  The  effect  of  the 
angle  0  between  the  ray  and  the  magnetic  field  was  very  small  and,  consequently, 
only  a  single  curve  has  been  used  to  represent  an  angle  of  10°  or  30°.  The  virtual 


Figure  4.  Variation  of  Virtual  Height  h'(0  With  Frequency  for  the 
Two  Parabolic  Models  (the  virtual  height  has  been  normalized  to 
maximum  real  height) 
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Figure  5.  Comparison  of  Exact  and  Calculated  Values 
of  Real  Height  for  the  Two  Parabolic  Models  (the  real 
height  has  been  normalized  to  maximum  real  height) 


height  fi.r  .Model  2  initially  rises  more  rapidly  than  that  for  Model  1,  then  lags 
behind  until  near  the  maximum  frequency  where  it  increases  very  sharply.  This 
behavior  in  virtual  height  may  bo  anticipated  qualitatively  from  the  electron  density 
as  a  funetmn  of  height  (see  solid  curves  in  I-'igure  a). 

A  comparison  of  the  exact  and  calculated  values  of  real  height  as  a  function  of 
input  frequency  is  shown  in  Figure  5.  In  view  of  the  several  approximations 
utilized  in  the  method,  the  agreement  is  generally  quite  good.  Model  1  showing 
slightly  better  results.  It  is  felt  that  the  less  accurate  results  for  Model  2  are 
due  primarily  to  numerical  inaccuracies  in  the  process  of  integration.  In  addition, 
the  occurrence  of  a  maximum  electron  density  for  Model  2  renders  the  calculated 
values  of  real  height  very  sensitive  to  small  errors  in  the  numerical  procedure 
near  the  maximum  height. 

As'  indicated  in  Section  5,  the  accuracy  of  the  method  may  be  improved  by  the 
introduction  of  additional  parameters.  However,  the  gain  in  accuracy  will  be  small 
while  the  additional  amount  of  computational  labor  required  will  be  quite  large. 

The  present  procedure,  utilizing  only  one  parameter,  offers  a  fairly  simple  method, 
without  the  burden  of  a  great  deal  of  numerical  work,  for  determining  real  height 
from  virtual  height  data. 

7.  SUMMARY  AND  CONCLUSIONS 

A  method  previously  developed  for  the  determination  of  real  height  from  virtual 
height  data  has  been  modified  and  simplified  by  representing  the  entire  dispersion 
curve  by  a  single  power  law  curve  with  only  one  parameter.  The  resulting  real 
height-virtual  height  relation  can  now  be  solved  in  a  simple  way  for  the  real  height. 
To  test  the  accuracy  of  the  method,  virtual  heights  were  computed  for  two  para¬ 
bolic  models  of  the  ionosphere  and  the  corresponding  real  height  values  then  ob¬ 
tained  by  the  present  method.  The  results  obtained  show  that  the  calculated  values 
of  real  height  were  in  very  good  agreement  with  the  exact  values. 

To  show  the  improvement  obtainable  with  the  introduction  of  additional  param¬ 
eters,  a  two-parameter  function,  which  gives  a  somewhat  more  accurate  represen¬ 
tation  of  the  dispersion  curve  then  does  the  one-parameter  function,  has  been  ob¬ 
tained  for  one  case.  However,  the  amount  of  computational  labor  required  is  con¬ 
siderably  greater  than  with  the  present  method. 

In  view  of  its  simplicity  and  accuracy,  it  is  felt  that  the  present  one-parameter 
method  is  a  useful  analytic  tool  for  determining  real  height  from  ionosonde  data. 
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